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ABSTRACT 


* 


The  Timation  n  satellite  Is  an  experimental  model  of  a 
range  and  doppler  measurement  satellite  navigation  system 
sponsored  by  the  Naval  Air  Systems  Command.  In  this  re¬ 
port  the  satellite  is  described  briefly,  along  with  the  ranging 
concept  and  the  range  and  doppler  measuring  equipment. 
The  basic  navigation  equations  for  range,  doppler,  and 
simultaneous  range  and  doppler  measurements  are  devel¬ 
oped,  along  with  several  alternate  techniques  depending  upon 
the  user  equipment.  An  analytical  expression  is  developed 
for  experimentally  determining  the  range  error  due  to 
first-order  ionospheric  refraction.  The  system  biases  are 
identified  with  clock  difference  (range  bias)  and  frequency 
difference  (doppler  bias).  An  error  analysis  is  performed 
to  determine  the  influence  of  geometry  and  system  noise 
upon  navigation  accuracy. 
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PRINCIPLES  AND  TECHNIQUES  OF  SATELLITE  NAVIGATION 
USING  THE  TIMATION  II  SATELLITE 


INTRODUCTION 

The  Timation*  experiment  for  satellite  navigation  is  being  developed  under  the 
sponsorship  of  the  Naval  Air  Systems  Command.  The  Timation  n  satellite  was  launched 
Sept.  30,  1969,  and  is  still  operating  at  this  date.  Timation  n  transmits  range  and  dop- 
pler  information  near  150  and  400  MHz,  which  can  be  used  to  correct  range  or  doppler 
for  first-order  ionospheric  refraction.  Five  U.S.-based  ground  stations  are  used  to 
track  the  satellite  and  collect  telemetry  information  from  the  sensors  on  board  the  sat¬ 
ellite.  Other  ground  stations  are  used  to  control  the  satellite' s  subsystems,  including 
the  ability  to  tune  (in  phase  and  frequency)  the  on  board  quartz  crystal  oscillator. 

The  overall  physical  configuration  for  Timation  n  is  given  in  Fig.  1.  Timation  n  is 
equipped  with  a  high-precision  quartz  crystal  oscillator  capable  of  frequency  stabilities 
on  the  order  of  a  few  parts  in  10 11  per  day.  Timation  n  is  equipped  with  active  thermal 
control  of  the  oscillator  environment,  which  effectively  eliminates  oscillator  frequency 
fluctuation  due  to  temperature  changes. 


Ranging  information  is  provided  by 
means  of  coherent  modulation  of  the  car¬ 
rier,  with  modulation  frequencies  varying 
from  100  Hz  to  1  MHz.  The  range  re¬ 
ceiver  synthesizes  a  similar  set  of  fre¬ 
quencies  which  are  phase  compared  with 
the  received  signal. 

The  doppler  measurements  are  ob¬ 
tained  by  adding  a  phase -locked  tracking 
filter  to  the  Timation  range  receiver. 
The  filter  is  used  to  track  the  doppler 
signal;  the  period  of  heterodyned  doppler 
is  averaged  over  a  specified  number  of 
cycles,  and  the  results  are  printed  for 
subsequent  date  analysis. 

In  this  report  the  equations  used  in 
the  analysis  are  developed,  and  the  biases 
of  the  range  and  doppler  measuring  sys¬ 
tems  are  Identified  with  the  time  differ¬ 
ence  (range  bias)  and  the  frequency  differ¬ 
ence  (doppler  bias)  between  the  satellite’ s 
and  observer's  clock  or  frequency  stand¬ 
ard,  A  sensitivity  analysis  is  performed 
to  determine  the  influence  of  geometry 
and  noise  on  navigation  accuracy. 


Fig.  1  -  The  Timation  satellite. 


"The  name  Timation  is  an  acronym  of  Time  Navigation. 
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THE  TIMATION  H  SATELLITE 

The  Timation  II  (Fig.  1)  satellite  has  an  overall  configuration  similar  to  the  Tima- 
tion  I  (1)  satellite;  hence  only  a  summary  of  the  features  will  be  given  in  this  report. 
The  satellite  weighs  approximately  125  lb  and  consumes  an  average  of  18  W  of  power 
furnished  by  solar  cells  and  batteries.  Two-axis  gravity  gradient  stabilization  is  pro¬ 
vided  by  using  an  extendable  boom.  Temperature  control  is  achieved  by  (a)  careful  de¬ 
sign  of  the  satellite  (2)  to  provide  a  temperature  range  from  0°C  to  +20° C  inside  the 
satellite  and  (b)  active  temperature  control  of  the  quartz  crystal  frequency  standard  to 
maintain  its  external  temperature  to  within  a  fraction  of  a  degree.  Linearly  polarized 
dipoles  are  used  for  the  150-  and  400-MHz  antennas.  A  separate  telemetry  antenna  is 
used.  This  is  mounted  on  the  side  and  has  more  than  40  dB  of  isolation  from  the  main 
antennas.  In  addition,  a  magnetometer  s  used  to  sense  attitude  changes  of  the  satellite, 

The  frequency  of  the  oscillator  may  be  electromechanically  tuned  in  discrete  steps 
of  approximately  3.6x10' 12 /pulse.  In  addition,  the  phase  of  its  transmissions  may  be 
advanced  or  retarded  in  discrete  steps  of  33.3  ns/pulse.  These  two  features  provide 
precise  control  over  the  satellite  clock  synchronization  and  clock  rate. 

Timation  n  is  in  a  500-naut-mi  near-circular  orbit  which  has  an  inclination  of  70° 
to  the  equatorial  plane.  With  this  orbit,  several  passes  of  12  to  16  min  duration  will  be 
available  during  the  day  at  each  of  the  five  Timation  tracking  stations. 


TIMATION  RANGING  CONCEPT 

The  Timation  n  satellite  carries  a  highly  stable  crystal  oscillator  from  which  nine 
modulation  frequencies  of  the  two  carriers  of  150  and  400  MHz  are  obtained.  The  modu¬ 
lation  frequencies  are  100  Hz,  312.5  Hz,  1  kHz,  3.125  kHz,  10  kHz,  31.250  kHz,  100  kHz, 
312.5  kHz,  and  1  MHz.  The  transmitted  modulation  frequencies  can  be  received  and 
phase  compared  with  a  similar  set  of  coherent  tones  synthesized  from  an  oscillator  or 
"clock"  at  the  receiver  site.  This  system  is  thus  a  frequency  interferometer  which  will 
measure  the  "time"  difference  between  the  received  signal  and  the  local  "time"  with 
ambiguities  of  60  ms  based  on  the  highest  common  divisor  of  12.5  Hz  and  an  accuracy 
based  on  the  precision  of  the  phase  comparison  of  the  highest  tone  (1  MHz).  In  the  sys¬ 
tem  the  resolution  of  the  phase  comparison  is  1%  of  a  period,  giving  a  time  resolution  of 
10  ns  when  using  the  1-MHz  "tone."  The  accuracy  of  the  "time"  comparison  of  the  re¬ 
ceived  and  local  signals  is  slightly  higher  than  10  ns  due  to  phasemeter  "zero"  adjust¬ 
ment,  nonlinearity,  differential  phase  shift  in  the  receiver,  noise,  etc.  This  measure¬ 
ment  may  be  converted  to  ranging  information  by  multiplying  by  c,  the  speed  of  light  in  a 
vacuum.  This  ranging  information,  which  depends  on  the  navigator's  position,  also  in¬ 
cludes  information  on  the  time  difference  between  the  satellite  clock  and  the  navigator's 
clock. 

The  actual  "time"  difference  between  the  received  signal  and  the  local  reference  is 
the  "time"  difference  between  the  satellite  oscillator  or  "clock"  and  the  ground  clock, 
plus  the  propagation  time  required  for  the  signal  to  propagate  from  the  satellite  to  the 
receiver.  The  "time"  indicated  by  the  components  of  the  received  signal  is  subject  to 
some  error  due  to  the  dispersive  effect  of  the  ionosphere. 

The  user's  time  base  is  obtained  from  the  user's  frequency  standard  using  suitable 
countdown  and  comparison  circuitry.  (The  results  of  time  transfer  via  Timation  n  will 
be  discussed  in  a  later  report.)  The  timing  requirements  for  the  ground  station  clock 
are  higher  than  for  the  user's  clock.  The  five  ground  stations  are  equipped  with  cesium- 
beam  frequency  standards  which  are  kept  in  time  synchronization  to  the  UTC  time  base. 
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Time  transfer  or  time  comparison  between  two  ground  clocks  can  be  made  by  com¬ 
paring  each  to  the  satellite  clock  and  combining  the  results.  There  are  two  conditions 
under  which  this  can  be  done:  (a)  the  case  in  which  locations  of  the  ground  clocks  are 
accurately  known,  and  (b)  the  case  with  unknown  locations.  These  two  cases  will  be  dis¬ 
cussed  later  in  the  report. 

The  system  user,  or  navigator,  is  not  required  to  have  the  same  precision  frequency 
standard  as  required  for  use  in  the  satellite.  For  example,  quartz  crystal  frequency 
standards  with  stabilities  on  the  order  of  a  few  parts  in  1010/day  would  be  suitable  for 
use  by  a  Timation  n  user.  Several  American  manufacturers  have  relatively-lcw-priced 
quartz  crystal  frequency  standards  which  meet  this  requirement.* 

A  practical  question  of  interest  is  how  a  navigator  would  initially  set  his  clock  and 
proceed  to  use  the  navigation  satellite.  A  simple  line  of  reasoning  will  be  used  to  show 
how  this  question  can  be  resolved. 

Figure  2  depicts  the  500-naut-mi  near- circular  orbit  obtained  for  Timation  n.  Ref¬ 
erence  to  Fig.  2  shows  that,  due  to  the  orbit  chosen  for  Timation  n,  the  range  for  a  navi¬ 
gator  on  or  near  the  Earth  must  be  between  500  and  1900  naut  mi,  which  corresponds  to 
propagation  times  varying  from  3  to  12  ms.  Now  for  the  purposes  of  discussion,  say  that 
the  satellite  transmits  ephemeris  information  on  the  minute,  which  includes  the  date, 
hour,  and  minute  of  transmission.  Two  cases  will  be  considered. 


Fig.  2  -  Schematic  of  500-naut-mi  near-circular 
orbit  obtained  for  Timation  II.  The  minimum  and 
maximum  ranges  from  an  observer,  on  or  near 
the  earth,  to  the  satellite  is  shown. 


Case  1  —  Navigator's  Position  on  Earth  Unknown 

Assume  that  the  navigator  is  receiving  information  from  the  satellite.  After  the 
satellite  transmits  the  first  minute  mark,  the  navigator  adjusts  his  clock  to  read  that 
minute,  plus  7.5  ms  (7.5  ms  is  the  average  transmission  time).  His  clock  will  now  be 
synchronized  to  within  ±5  ms  of  the  satellite  clock.  The  navigator  now  proceeds  to  take 


*See  Ref.  3  for  a  survey  on  the  recent  advances  in  frequency  standards. 
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data  until  the  satellite  goes  below  his  horizon.  He  should  now  note  the  measured  range 
for  the  last  minute  of  observable  data.  That  range  must  be  (see  Fig.  2)  near  12  ms; 
hence  another  correction  may  be  made  to  bring  the  navigator’ s  clock  to  within  ±1  ms  of 
the  satellite  clock.  Note  that  so  far  no  knowledge  of  the  navigator’s  position  has  been 
assumed,  except  that  he  is  on  or  near  the  earth. 

Case  2  —  Position  Approximately  Known 

Since  the  ranging  system  is  modulo  80  ms  and  the  time  is  known  to  ±1  ms,  the  rang¬ 
ing  system  may  now  be  used  to  provide  a  much  higher  precision  time  synchronization. 
This  is  done  by  navigating  while  using  all  of  the  ranging  information  obtained  as  the  sat¬ 
ellite  sweeps  overhead.  However,  to  do  this,  an  approximate  position  must  be  assumed 
to  solve  the  nonlinear  equations  involved.  The  exact  radius  of  convergence  required  for 
these  equations  is  not  known,  except  that  it  is  large  compared  with  any  reasonable  un¬ 
certainty  in  position.  For  example,  convergence  has  been  obtained  using  an  assumed 
position  of  1000  naut  mi  further  away  from  the  known  position.  Needless  to  say,  this 
should  be  sufficient  to  cover  most  situations  of  interest. 


TIMATION  n  RANGE  AND  DOPPLER  MEASUREMENT 

Figure  3  shows  a  simplified  block  diagram  of  the  Timation  n  400-MHz  range  re¬ 
ceiver.  The  Timation  II  satellite  is  tracked  using  a  high-gain  directional  antenna  fol¬ 
lowed  by  a  preamplifier  and  two  stages  of  IF  at  frequencies  of  370  and  25  MHz.  The 
received  tones  (varying  from  100  Hz  to  1  MHz)  are  then  detected,  using  a  double  detec¬ 
tion  technique  for  the  lower  tones,  and  phase  compared  with  a  similar  set  of  tones 
(synthesized  from  the  5-MHz  frequency  standard)  to  produce  the  phase  difference  which 
may  be  combined  to  yield  the  range  measurement. 
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Fig.  3  -  Block  diagram  of  the  Timation  II  400-MHz  range  receiver. 


Doppler  measurements  are  obtained  only  at  the  NRL  ground  station.  The  NRL 
ground  station  is  equipped  with  a  phase-locked  tracking  filter  to  measure  the  satellite’s 
doppler  shift  at  399.4  MHz.  This  is  accomplished  by  heterodyning  a  4.38-MHz  signal 
(synthesized  from  the  5-MHz  standard)  with  the  4.4-MHz  signal  available  as  a  result  of 
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the  two  stages  of  mixing  used  to  derive  the  ranging  information.  The  procedure  results 
in  the  doppler  signal  being  offset  by  a  constant  20  kHz  The  phase-locked  loop  is  then 
used  to  track  the  signal,  which  is  counted  and  printed  for  subsequent  data  reduction. 


SATELLITE  TRAJECTORY  CALCULATIONS 

The  satellite  trajectory  computation  is  made  by  the  Naval  Weapons  Laboratory 
(NWL)  using  doppler  tracking  data  obtained  from  the  NWL  tracking  network  (TRANET). 
Although  the  range  data  from  the  five  Timation  ground  stations  has  been  included  in  the 
orbit  determination,  it  is  not  routinely  included  because  of  the  data  communications  re¬ 
quired  of  the  NWL  operating  system. 

The  orbit  determination  is  performed  on  the  NWL  STRETCH  computer  using  their 
ASTRO  (4)  program  which  performs  a  statistical  estimate  of  the  dynamical  and  observa¬ 
tional  parameters  of  the  state  variables  at  epoch.  The  force  model  accounts  for  accel¬ 
erations  from  the  following  sources:  (a)  earth  gravitational  accelerations,  (b)  sun  and 
moon  gravitational  accelerations,  (c)  solar  and  lunar  tidal  bulge  effects,  (d)  atmospheric 
drag,  and  (e)  radiation  pressure.  The  earth’s  gravitational  acceleration  includes  coeffi¬ 
cients  for  the  earth's  gravitational  potential  as  a  function  of  longitude  as  well  as  latitude. 
Other  parameters  such  as  drag  and  the  position  of  the  tracking  stations  are  included  in 
the  model.  A  weighted  least-squares  estimate  is  then  performed  based  on  observational 
data  obtained  over  time  arcs  ranging  from  2  to  5  days. 

The  Timation  II  satellite  is  equipped  with  dual  frequencies,  which  allows  lor  the 
measurement  of  first-order  ionospheric  refraction.  With  the  inclusion  of  the  ionospheric 
refraction,  NWL  determines  the  position  of  Timation  II  to  ±20  m  during  the  observation 
span.  The  positional  accuracy  outside  of  the  observed  data  span  remains  near  ±20  for 
extrapolations  on  the  order  of  12  to  24  hr.  Beyond  1-day  extrapolations,  the  error  may 
grow  rapidly. 

For  operational  purposes  the  satellite  ephemeris  would  require  updating  on  a  fre¬ 
quent  basis.  However,  for  the  purpose  of  analyzing  the  Timation  system  performance, 
the  analysis  is  routinely  done  using  the  satellite  trajectory  during  the  observed  data 
span.  This  choice  minimizes  the  contribution  of  satellite  positional  error  to  the  overall 
system  error. 


COORDINATE  SYSTEM 

The  satellite  trajectory  is  given  with  respect  to  an  earth-fixed  coordinate  system. 
Figure  4  depicts  the  coordinates  of  the  satellite  and  the  observer  with  respect  to  this 
frame.  The  z  axis  is  the  earth’s  polar  axis,  the  x  axis  is  along  the  Greenwich  meridian 
in  the  equatorial  plane,  and  the  y  axis  completes  the  right-handed  set.  The  satellite's 
position  r  and  velocity  v  with  respect  to  this  set  is  given  by 

r  =  (x,y,z) 

and 


V  =  (Vx'Vvz>  ■ 

The  observer's  coordinates  with  respect  to  this  set  is  denoted  by  r8  =  (xs.ys.zs).  The 
slant  range  between  the  observer  and  the  satellite  is  given  by 


R  =  lr-  rsl 


(1) 
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Fig.  4  -  Earth-fixed  coordinate  system  used 
to  locate  the  satellite's  position  r  and  the 
observer  r,. 


The  rate  of  change  of  R  is  denoted  as  R ,  and  a  unit  vector  along  the  slant  range  is 
given  by 


The  satellite  ephemeris  time  base  will  be  denoted  by  ts  and  the  observer's  time 
base  is  denoted  by  t  with  the  two  bases  related  by 

ts  =  t,  +  At  .  (3) 

The  quantity  At  should  not  be  confused  with  the  transmission  time  of  the  RF  signal 
from  the  satellite  to  the  observer.  The  time  transmission  delay  is  given  by  -R/c,  where 
e  is  the  speed  of  light  in  vacuo.  A  correction  for  the  time  transmission  delay  may  be 
applied  to  either  the  measured  data  or  the  ephemeris.  An  explicit  term  for  this  will  not 
be  shown  in  the  equations  developed  in  the  text,  but  the  reader  should  keep  in  mind  that 
this  correction  should  be  accounted  for. 


RANGE  NAVIGATION  EQUATIONS 

The  basic  range  navigation  equations  reflect  the  three  quantities  (a)  time  difference 
At  between  the  observer’s  clock  and  the  satellite's  clock,  (b)  the  observer's  latitude, 
and  (c)  the  observer's  longitude.  The  frequency  difference  between  the  satellite  and 
observer  frequency  standards  is  measured  using  values  of  At  obtained  from  successive 
satellite  passes  taken  from  different  positions.  Denote  the  measured  range  by  R0  (the 
system  measures  time  difference,  which  is  readily  converted  to  range  by  multiplication 
by  c ,  the  velocity  of  light)  and  the  computed  range  (using  an  assumed  position  rs  for  the 
observer)  by  Rc.  Then  the  observed  minus  the  computed  time  difference  for  each  point 
measured  during  the  satellite  pass,  denoted  by  (0-  C)j,  is  given  by 

(O-C)  (R0)  -  (K  +  (Rc)  )  .  (4) 

i  i  i 

In  Eq.  (4),  K  c  (At )  where  At  ts  -  t  is  the  actual  time  difference  between  the  satel¬ 
lite's  and  the  observer's  clocks.  Now,  a  least- squares  fit  is  preformed  which  adjusts 
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these  three  quantities  in  such  a  way  that  the  sum  of  squares  of  the  range  residuals  are 
minimized. 

The  explicit  function  to  be  minimized,  denoted  by  <£( K,  x,  #),  is  given 


*(K.  K,9)  =  £  (°-c)i2 


S  {(ba  -  [ 


(Rc).  +  K 


(5) 


In  Eq.  (5),  n  is  the  number  of  range  observations  made  during  a  satellite  pass,  and 
the  computed  function  is  =  K  +  (Rc)  . 

The  technique  used  to  solve  Eq.  (5)  for  a  minimum  is  the  least-square  procedure  (5) 
which  involves  the  linearization  of  the  computed  function  about  assumed  values  of  K,  X, 
and  x.  In  matrix  form,  the  least-square  solution  is  given  by 

BAP  =  E  (6) 


where 


E  =  ATW  [0  -  C  ]  ,  and  B  r  ATHA  . 

The  3x1  column  vector  ap  has  components  Ak,  ax,  a e  which  are  the  corrections  to 
the  assumed  values  of  K,  X,  d,  the  N  x  1  column  vector  o-C  has  the  range  residuals 
( R0 ) .  -  [(  Rc ) .  +  k]  as  components,  A  is  the  N  x3  matrix  of  partials  of  the  computed  func¬ 
tion  Cj  (with  respect  to  K,  X,  e)  evaluated  for  each  observation,  w  is  the  nxn  weight 
matrix,  and  at  is  the  a  matrix  transposed. 

The  required  partial  derivatives  for  the  matrix  A  are  3C;/bk,  BCj/BX,  and  BCj/B#. 
The  observer's  position  rs  is  related  to  the  latitude  and  longitude  by  the  spherical-to- 
rectangular  coordinate  equations  given  by 

xs  =  ^a  cos^  COsX 

ys  =  Rs  cos 6  sinX  (7) 

z  =  R  sin#  . 
s  s 

where  e  is  the  geocentric  latitude  and  X  is  the  longitude  measured  East  from  the  Green¬ 
wich  meridian. 

The  required  partials  are  given  by 


BCi 

3K 

'  1 

3C. 

3(RC). 

I 

Bxs 

B(Rc). 

1 

--R  -  ^ 

BX 

Bxs 

~3X 

BX 

BX 

3C. 

3  (RC  )  i 

dxa 

‘d(Rc)i 

Bys 

+  3{Rc)i  Bzs 

b  e 

9x. 

b  e 

3vs 

3# 

+  dz,  d 6 

3# 


(8) 
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Fig.  5  -  In  one  formulation  of  the  range  navigation  equations,  a  constraint 
c  is  defined  (see  Eq.  13).  This  figure  shows  that  the  gradient  of  c  is  a 
vector  perpendicular  to  the  surface  of  the  earth. 


Figure  5  illustrates  that  the  equation  of  condition  is  equivalent  to  saying  that  the  dot 
product  of  rfl  with  a  vector  7g  (gradient  of  g)  normal  to  the  surface  at  the  assumed  posi 
tion  of  the  observer  equals  zero. 

A  third  useful  formulation  is  obtained  by  constraining  the  solution  to  the  plane  deter 
mined  by  the  velocity  vector  at  the  time  of  closest  approach  (denoted  by  TCA)  of  the 
satellite  and  the  radius  vector  from  the  observer  to  the  satellite  at  TCA. 


COMPONENT  OF  Ar,  IN 
PLANE  DEFINED  8Y 


Fig.  6  -  Arrangement  of  the  pertinent  range  navigation  vectors  when  the  range  solution 
is  confined  to  the  plane  defined  by  »TCA  and  vTcA. 


Figure  6  shows  that  this  procedure  is  equivalent  to  decomposing  A<?,  A\  into  two 
components,  one  in  the  observer's  horizontal  plane,  and  the  other  a  projection  onto  the 
radius  vector  at  the  TCA. 
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A  vector  perpendicular  to  the  plane  defined  by  VTCA  and  (r  -  r,)TC*  is  given  by 
T  -  (r  -  rs)TCA  x  VTC*.  This  vector  is  constant,  using  the  assumed  value  of  r„;  for  an 
arbitrary  rs  the  constraint  is  expressed  by 


h  =  r.  T  (15) 

The  partials  bh/Bx,, ,  BH/3yg,  and  3H/3Zg  are  required  in  the  least-square  formulation 
and  are  given  by 


BH 

3r. 


BH 

BH 

*y. 

bh 


(y-  y.>TCAVzTCA  -  "  **)TCAVyTCA 

"  (  *  ~  xi  )TCA  VzTCA  *  (  *  “  1  ■  >TCA  VxTCA 
(x_  x«)'fCAVyTCA  *  <y  y.>TCA  VxTCA 


(16) 


The  matrix  shown  in  Eq.  (17)  below  gives  the  solution  which  will  be  referred  to  as 
the  "along  track"  solution.  Using  VH  obtained  from  Eq.  (16),  the  last  equation  in  the 
matrix  of  Eq.  (17)  expresses  this  constraint: 


(17) 


The  corrections  Ax  ,  Ayt,  Azs  may  be  decomposed  into  components  along  the  track 
of  the  satellite  (ARV)  ana  along  the  vector  Rtca  which  is  denoted  by  ARIt.  This  decompo¬ 
sition  is  useful  for  studying  the  effects  of  several  error  sources  on  the  navigation  solu¬ 
tion,  some  which  will  be  discussed  later.  The  final  latitude  and  longitude  may  be  obtained 
by  correcting  x„  y„  z,  and  converting  from  rectangular  to  spherical  coordinates,  or  by 
obtaining  the  two  components  in  the  horizontal  plane  and  then  applying  the  parallelogram 
law  to  obtain  the  (vector)  correction  to  the  observer's  assumed  position. 
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The  magnitude  of  ARR  and  ARV  is  given  by 

|ARV|  ,  •  Ar.  (18) 

' VTCA I 

and 

|  ARR|  kfCh  •  Ar,  (19) 

The  three  formulations  listed  above  do  not  exhaust  the  possible  ways  to  obtain  a 
position  fix  from  the  range  data.  Notice  that  in  Eq.  (4)  the  K  is  assumed  constant  during 
the  pass  (up  to  16  min  duration).  The  reason  for  this  is  that  the  Timation  System  is  de¬ 
signed  to  utilize  the  high-precision  (3,6)  quartz  crystal  frequency  standards  available 
today.  It  should  be  emphasized  that  it  is  the  uncertainty  in  the  frequency  that  is  impor¬ 
tant,  inasmuch  as  the  frequency  offset  due  to  crystal  aging  can  be  computed  as  a  known 
correction  to  the  measured  data. 

It  is  not  necessary  that  the  observer's  frequency  be  precisely  known  in  order  to  get 
a  fix  using  Timation  n.  A  different  computed  function  may  be  used  where  the  observed 
data  is  the  range  difference  between  two  successive  minutes.  Note  that  continuous 
measurement  of  the  range  between  minutes  is  not  required  to  obtain  this  data.  The  navi¬ 
gation  equations  needed  to  utilize  data  in  this  fashion  are  related  to  navigation  equations 
obtained  by  integrating  the  expression  for  doppler,  which  is  sometimes  called  the  "inte¬ 
grated  doppler"  solution.  The  relationship  between  "range  difference"  and  "integrated 
doppler"  will  be  discussed  after  the  doppler  navigation  equations  are  derived. 

What  will  be  shown  is  that,  for  data  treated  in  this  manner,  the  clock  difference 
(which  may  be  expressed  in  terms  of  range  by  multiplication  by  c,  the  speed  of  light)  is 
eliminated  from  the  observed  range  data,  and  the  bias  parameter  is  the  product  of  the 
navigator' 8  frequency  difference  by  the  difference  in  observation  time.  The  navigator's 
frequency  difference  is  related  as  a  constant  rate  of  change  of  range  through 

",  =  -(f;)(f‘- f«)  •  (20) 

Let  the  measured  ranges  for  each  minute  be  denoted  by  (R,m,  t, ),  (R™,  t,) . . .  (RNm,  tN) 
where  ( t  a  -  t  , )  ■  at;  ( t ,  -  t , )  ■  AT ,  etc.  Let  rb  be  the  rate  of  change  of  range  due  to  the 
difference  ( fs  -  fg)  between  the  satellite  frequency  standard  and  the  navigator's  fre¬ 
quency  standard.  Let  R_,  represent  the  range  difference  due  to  the  difference  At  - 
t,  -  t,  between  the  satellite  clock  and  the  navigator's  clock  at  time  t, .  Then  the  meas¬ 
ured  ranges  R^  may  be  separated  to  isolate  the  parts  due  to  At  and  ft  -  fg  by 


K  R,  *  R«I  ♦  K  <*I-  V 

Rjm  -  RJ  '  R*I  +  Rg  («a- 

K  1  R3  ♦  R,i  ♦  V'j*  *«> 

RNm  -  RN  +  Rgl  '  V'N-  »1> 


(21) 
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The  measured  ranges  may  be  differenced  to  obtain 

R2m  -  R™  =  (R2  -  R, )  +  RgAT 

Rjm  -  R2m  =  (R3-  R2)  +  RgAT 

R4m  -  R™  =  (R4-R3)  +  RgAT  (22) 


RfT  “  RN-  I  =  ~  RN-  i  )  +  Rg  ■ 

Note  that  by  differencing,  the  term  Rgl  which  represents  the  clock  difference  At  (ex¬ 
pressed  in  terms  of  range)  has  been  eliminated.  A  new  computed  function  may  now  be 
defined  as  c  =  k  +  |R(  t  r  at)  -  r<  t )  ] ,  where  k  =  Rg  at.  The  navigation  equations  for  the 
"range  difference"  will  be  discussed  after  development  of  the  doppler  navigation  equa¬ 
tions. 

The  previous  development  illustrates  that  the  range  data  may  be  processed  with  the 
frequency  difference  (f  -  fg)  known  or  unknown.  The  decision  as  to  which  model  the 
navigator  should  use  will  be  influenced  by  the  stability  of  his  clock.  It  should  be  empha¬ 
sized  that  if  a  user  navigates  using  k  =  ia  at,  the  clock  difference  At  =  t,  -  tg  may  be 
recovered  after  the  user's  position  has  been  determined  by  evaluating  R™  -  R,  +  rrI  , 
with  At  =  Rfl/c. 


DOPPLER  NAVIGATION  EQUATIONS 

The  basic  doppler  navigation  equations  reflect  the  three  quantities  (a)  frequency 
difference  between  the  observer's  frequency  standard  and  the  satellite  frequency  stand¬ 
ard,  (b)  the  observer's  latitude,  and  (c)  the  observer's  longitude.  The  fourth  required 
quantity,  At,  the  difference  between  the  observer's  clock  and  the  satellite  ephemeris, 
is  assumed  to  be  known.  The  least-squares  technique  will  again  be  used  to  obtain  a 
solution;  however  the  function  to  be  minimized  is  the  sum  of  the  squares  of  the  frequency 
residuals. 

The  computed  doppler  shift,  denoted  by  dopcJ  ,  is  given  by  Eq.  (23),  where  f„  de¬ 
notes  the  satellite  system  frequency. 


(23) 


Equation  (23)  may  be  rewritten  to  show  the  explicit  dependence  of  doppler  shift  upon  the 
position  of  the  observer: 


({.\ vxi(xi  - *,)  1  vyi(yj-y.)  + 

.  i  \  c  /  j  -  -  y 

V(Xi  -  X,)  +  (Yi  -  ys)  +  (*i  -  *,) 


.  Note  that  in  Eq.(24)  no  term  is  present  to  reflect  the  observer's  velocity  with  re¬ 
spect  to  the  earth-fixed  set.  This  could  be  accomplished  by  replacing  vc ,  by  vcl  -  vs  ( t ) 
where  vs(t)  is  the  observer's  velocity  (known  a  priori). 


I' 
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The  computed  function  is  given  by 

Cs  K  *  DOPci  (25) 

where  K  now  represents  the  frequency  difference  (K  assumed  constant). 

The  observed  minus  computed  function  is  given  by 

(O  -  C).  DOPoi  -  (K  t  DOPc  j )  (26) 

The  partial  derivatives  required  to  obtain  a  solution  are  lengthy,  hence  only  one 
formulation  will  be  given.  The  required  partials  for  this  formulation  are  given  by 


3DOPci 

G?) 

Vxl 

l  mp 

(Xi  -  X,) 

R1 

t  UOPci 

R>  ' 

aDOPci  = 

(M 

vyi 

4  nnp 

<y,  -  y,) 

Ri 

+  UUKci 

R  2 

K1 

9DOPci 

(M 

v,. 

i  rmp 

(Z1  "  z») 

Bzf 

\  C  / 

R! 

+  IAJKci 

RiJ 

The  matrix  equation  shown  below  gives  the  least-squares  solution,  which  is  a  con¬ 
strained  solution  for  the  latitude  and  longitude  as  a  function  of  xs,  yt ,  zt: 


(28) 


Doppler  solutions  in  terms  of  the  bias  parameter  K,  the  along-track  error  ARV,  and  the 
slant  range  error  ARR  are  obtained  in  addition  to  the  solution  for  latitude  and  longitude. 

The  measured  doppler  shift  dop0  is  obtained  from  the  1  k  (1000  cycles),  10  k,  or 
40  k  period  average  data  available  with  some  of  the  Timation  n  equipment.  To  obtain 
the  conversion,  let  t0  be  the  midpoint  of  the  measurement  interval,  N  the  number  of 
cycles  counted  by  the  equipment,  and  tn  the  total  amount  of  time  required  to  count  N 
cycles  of  the  offset  doppler  signal.  The  ground  station  time  base  tK  is  used  to  obtain 
the  count  time  tn  .  Let  f  represent  the  frequency  that  is  counted,  which  is  given  by 
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f  =  R(t0)  ♦  (q  *  (fn~  fK>  (29) 

where  the  first  term  is  the  doppler  shift  at  t0,  fQ  is  a  constant  offset  in  the  ground  sta¬ 
tion  equipment,  and  (a  -  f„  represents  the  difference  between  the  satellite  and  ground 
station  frequency  standards.  The  number  of  cycles  n  obtained  by  counting  f  over  a  time 
interval  tn  centered  about  point  1 0  is  given  by 


*o*T*/2  f 

N  -  J  fdt  =  [R(t0+  Tn/2)  -  R(t0-V2)|  ♦  K0  +  (fs  -  fR))TN 

«o-Tn/2 


(30) 


Since  R  may  be  expanded  about  the  point  t0,  Eq.  (30)  may  be  approximated  by 

N  =  (•■  'i  R  ♦  ( f<j  +  K)  +  R  ^  TNJ)  Tn  •  (3D 

The  term  R  (tQ>  tn2  /24  is  the  measure  of  the  error  which  is  committed  in  obtaining 
the  measured  doppler  shift  dop0  given  by 


DOP,  .(£)-  f0  (32) 

\ 

Figure  7  is  a  plot  of  the  maximum  error  committed  by  using  Eq.  (32)  to  obtain  the 
measured  doppler.  However,  a  correction  has  been  applied  in  the  results  for  this  term 
so  that,  for  the  longest  count  used  40,000  cycles,  the  representation  error  will  be  con¬ 
siderably  smaller  than  the  measurement  error. 


Fig.  7  -  Maximum  error  committed  by  using  the  relation 
DOP0  =  (N/T)  -  fQ  for  the  measured  doppler. 
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A  navigation  may  also  be  performed  using  Eq.  (30),  which  is  sometimes  called  the 
"integrated  doppler"  solution.  Comparison  of  Eq.  (30)  with  the  range  difference  Eq.  (21) 
shows  .the  relationship  between  the  "integrated  doppler"  and  the  range  differences  ob¬ 
tained  by  differencing  the  range  data.  The  factor  -  f  /c  is  a  constant  and  the  extra  term 
f  p  tn  in  Eq.  (30)  is  constant  whenever  tn  is  specified  by  the  user.  The  factor  fQ  is  de¬ 
pendent  on  the  counting  equipment  used,  and  the  NHL  ground  station  uses  an  offset  of  20 
kHz.  The  term  ( f  s  -  fg )  tn  may  now  be  related  to  the  term  Rg  in  the  range  difference 
solution  by 


It  can  now  be  seen  that  Eqs.  (30)  and  (21)  are  related  even  though  the  measurements 
of  "range  difference"  and  "integrated  doppler"  are  obtained  by  two  independent  measur¬ 
ing  teclmiques.  As  stated  previously,  the  quantity  At  -  t5  -  tg  cannot  be  obtained  by  the 
"integrated  doppler"  technique,  and  therefore  the  clock  correction  cannot  be  obtained  by 
this  method. 


The  navigation  equations  for  "integrated  doppler"  will  be  derived,  with  the  "range 
difference"  equations  obtained  by  using  Eq.  (30),  multiplying  the  first  term  by  -c/ft , 
and  removing  the  term  fQTN  from  the  computed  function. 

The  auxiliary  function  D  ( t ,  tn  ) ,  given  by 

iR(t„  +  TN/2)  -  R(t0- yaji  , 


is  defined  for  convenience,  and  K  -  ( f s  -  fg)  TN  will  be  the  bias  parameter. 
Now  the  computed  function  is  given  by 


C.  =  K  + 

The  reqi  red  partials  are  given  by 


BC; 

3K 

=  i 

BC. 

1 

-L  = 

_  f. 

BR 

Bx 

9 

Bx 

s 

c 

9D; 

"  BR 

Bys 

c 

.3y" 

*2 

3D. 

i 

'  BR 

9zs 

c 

3z, 

D»(trTN^  f 


(ti  + Tn/2) 

(*i  +  V2) 
(t.  +Tn/2) 


fQTN 


<‘i-V2> 


(ti- V2> 


(t.  -Tn/2) 


(33) 


(34) 


The  equations  given  will  be  using  the  constraint  given  in  Eq.  (35)  below.  However, 
the  along -track  solution,  or  the  direct  solution  for  latitude  and  longitude,  may  be  obtained 
in  a  manner  similar  to  the  technique  used  for  the  corresponding  range  solutions. 


The  frequency  difference  f,  -  fR  may  be  determined  if  AT  is  a  variable,  that  is,  the 
product  ( f,  -  f_)  AT  is  not  constant.  This  may  be  done  by  setting  C|  =  (fs-  f,)ATi 
*  d  ( t . .  T.)  and  using  the  constant  term  as  f  s  -  fg.  This  may  be  done  by  replacing  ac.  / 
3k  l  by  SCj/Bcf,-  fg)  AT .  in  the  normal  equations. 
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This  concludes  the  development  of  the  navigation  equations  for  range,  range  differ¬ 
ence,  doppler,  and  integrated  doppler  measurements.  A  data  weighting  matrix  may  also 
be  used  according  to  the  techniques  given  by  Ref.  (5);  however  this  was  not  done  in  this 
development  in  order  to  simplify  notation.  Other  iterative  techniques  exist,  and  it  should 
be  mentioned  that  graphical  techniques  (7)  may  be  used  to  obtain  solutions  for  either 
range  or  doppler  position  fixes. 


RANGE/DOPPLER  NAVIGATION  EQUATIONS 

The  simultaneous  measurement  of  both  range  and  doppler  offers  the  potential  of  an 
observer  obtaining  a  two-dimensional  fix  (latitude  and  longitude)  "instantaneously."  In 
order  to  do  this  the  range  and  doppler  bias  parameters  must  be  known.  This  require¬ 
ment  indicates  than  an  atomic  frequency  standard  would  be  desirable  for  the  user.  The 
atomic  frequency  standards  commercially  available  are  capable  of  extremely  stable 
phase  and  frequency  operation.  Hence  a  system  user  equipped  with  an  atomic  frequency 
standard  could  place  a  close  tolerance  on  both  the  range  and  frequency  bias  parameters 
for  a  relatively  long  interval  of  time  after  some  initial  calibration.  Note  also  the  fact 
that  some  of  the  ultrahigh-precision  quartz  crystal  oscillators  may  also  be  used,  but  for 
shorter  intervals  of  time  between  calibrations,  which  may  be  accomplished  by  periodi¬ 
cally  taking  more  data  than  the  two  points  required  for  an  "instantaneous  fix." 

The  three  equations  that  are  used  to  obtain  a  solution,  using  an  assumed  position, 
are  given  by 


I  **  -  rg  |  =  R  , 

(r  -  r8) 


DOP  , 


r  -  r. 


and 


'.I  =  Rs 


where  the  first  two  equations  are  repeats  of  Eqs.  (1)  and  (23). 

These  three  equations  are  rewritten  slightly  for  convenience  as 
f s 

-~T  (v  (x-  x,)  +  v  (  y  -  y  )  +  v_  (  z  -  zj]  =  R(DOP) 


(36) 


(37) 


(X  -  Xs)2  +  (y  -  ys)2 


( z  - 


R2 


and 


+  y. 


R  2 

S 


(38) 

(39) 


The  observed  quantities  are  denoted  by  R0  and  DOP0 .  There  is  no  need  to  use  least 
squares  to  get  a  solution  here  because  we  have  exactly  three  equations  of  condition  in 
three  unknowns  (xs,  ys,  z5).  The  differential  correction  equations  are  given  by  the 
matrix 
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f  f  f 

r  i 

/f\ 

n  s  s 

—  v  —  v  —  v 

c  x  c  y  c  z 

DOP0Ro  -  l— J  v  •  (r-  rs) 

-(X  -  Xs)  -(y-  ys)  -(z  -  zs) 

AVs 

= 

(R02  -  R2X2 

xs  zs 

Azs 

0 

(40) 


It  should  be  noted  that  a  polynomial  solution  in  one  of  the  three  coordinates  can  be 
obtained  by  suitable  manipulation.  Undoubtedly,  a  trigonometric  solution  in  terms  of 
latitude  or  longitude  could  be  obtained  by  expressing  the  range  and  doppler  as  functions 
of  6  and  \.  However,  Eq.  (40)  is  relatively  easy  to  compute  and  converges  well  for  any 
reasonable  guess  of  the  observer's  position  rs. 


FIRST-ORDER  RANGE  IONOSPHERIC  CORRECTIONS 

The  ranging  information  at  150  and  400  MHz  may  be  combined  to  measure  and  re¬ 
move  first-order  ionospheric  refraction  effects.  In  Fig.  8  denote  the  actual  transmis¬ 
sion  times  for  the  150-  and  400-MHz  signals  by  T1S0  and  T400,  respectively.  Equations 
(41)  and  (42)  give  an  expression  for  the  total  path  length,  where  n  is  the  index  of  refrac¬ 
tion  (along  the  optical  path),  ds  is  the  element  of  arc  length,  and  pt  and  PR  denote  the 
position  of  the  satellite's  transmitting  antenna  and  the  observer's  receiving  antenna, 
respectively: 


T 


ISO 


(41) 


T 


400 


(42) 


Although  the  path  of  integration  is  slightly  different  for  each  signal,  each  path  i  s  close  to 
the  geometric  path  along  R . 


Fig.  8  -  Typical  ray  path  trajectory.  The  quantities  and  t400 
give  the  actual  transmission  times  for  the  150-  ana  400-MHz 
ranging  signals,  respectively. 
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The  index  of  refraction  will  be  split  into  two  parts,  with  An,  rop  representing  the 
error  (An,rop  n,  rop  -  l)  due  to  the  troposphere,  which  is  assumed  to  be  independent 
of  frequency.  That  is,  n'1S0  n,SQ  *  Antrop,  and  n'00  -■  n400  f  A„trop. 

The  ionospheric  refraction  may  be  expanded  in  the  form  given  by 

a  P 

n  ‘  1  “  Fa  +  71  ■  •  ■  (43) 


with  a  and  p  dependent  upon  time  and  position. 


Let  f ,  denote  the  150-MHz  frequency  and  i  2  denote  the  400-MHz  frequency.  Also 
note  that  Cf ,  =  3  f2.  The  actual  ratio  (3/8)  is  slightly  different  from  an  exact  ratio,  but 
the  form  of  the  correction  may  be  easily  changed  to  accommodate  any  ratio.  Substituting 
Eqs.  (41)  and  (42)  into  (43),  and  noting  that  for  n  =  1  we  obtain  the  geometric  path  from 
pt  to  PR ,  we  obtain 


and 


:  -  f 

150  J 


A,1tropris 


-  r  t?  *  r 

i 


Pds 

"v 


(44) 


T 


4  00 


rp* 

J  A"tropds 


R 

c 


(45) 


Now  using  the  relation  8  f ,  =  3  f2,  the  first-order  ionospheric  term  is  eliminated, 
resulting  in 


R 

c 


9_ 

55 


^4  0  0  _  5  0  ) 


rPR 

-  J  Antropfls 

PT 


p 

320  f  R  PA s 
165  J  -fj 

P  T  1 


(46) 


Let  AR1S0  and  ar400  denote  the  first-order  corrections  to  the  150-  and  400-MHz  ranges, 
respectivei> .  Then 


AR 


AR 


“oo  -  55 

64  c 


5  ^T400  ^150^ 


32c 


1  5  0 


~51>  ^T400  _  1  5  0  )  * 


8 2  -  32 
8  2c 

8 2  -  32 


(^400  _  50) 

(^4  00  '  "^150^ 


(47) 

(48) 


The  tropospheric  term  in  Eq.  (46)  can  be  calculated  using  measurements,  or  a 
standard  atmosphere  can  be  used  to  obtain  an  approximate  value. 

The  form  of  Eqs.  (47)  and  (48)  is  such  that  it  appears  to  be  easier  to  make  the  cor¬ 
rection  by  computer  rather  than  by  additional  circuitry. 


SINGLE -PASS  GEOMETRY 

An  important  consideration  to  a  system  user  is  the  geometry  presented  by  a  system 
configuration  and  the  influence  of  system  geometry  and  overall  system  noise  to  naviga¬ 
tional  accuracy.  The  range  and  doppler  navigation  equations  derived  in  this  report  may 
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be  used  for  data  gathered  in  one  pass  of  the  satellite,  which  may  be  as  much  as  16  min 
in  duration.  Reference  7  gives  other  possible  satellite  configurations  for  Timation. 

A  weighted  least-squares  analysis  may  be  made  with  the  data  which  includes  such 
factors  as  the  uncertainty  in  the  estimate  of  the  assumed  position  and  the  correlations 
present  in  the  observed  data  with  appropriate  weighting  factors.  Such  a  weighted  least- 
squares  formulation  may  be  obtained  using  techniques  derived  in  Ref.  8.  Equation  (49) 
below  gives  one  such  formulation  where  f  x  is  the  a  priori  state  covariance  matrix,  rx  is 
the  "post-flight"  covariance  matrix,  A  is  the  matrix  of  observational  partials  as  defined 
in  this  report,  w  is  a  weight  matrix,  and  rx  is  the  data  covariance  matrix: 

rx  =  (ATWA+ (ATwrf  WA  +  f-,)(ATWA+ fx)_1  .  (49) 

The  geometry  may  be  singled  out  for  study  by  assuming  that  (a)  rx_  1  =  o ,  (b)  rx  is 
uncorrelated  with  r;1  =  w,  and  (c)  w  =  (l /a2)  I.  With  these  assumptions,  Eq.  (48)  sim¬ 
plifies  to 


rx  r  CT2(  ATA)  ”  1  =  cr2B"1  .  (50) 

Now  Eq.  (50)  may  be  evaluated  for  a  selected  single-pass  geometry,  using  the  appropri¬ 
ate  equations  for  range  or  doppler.  Figure  9  indicates  that  for  the  single-pass  geometry 
the  important  factors  are  (a)  maximum  elevation  angle,  (b)  symmetrical  data,  and  (c) 
amount  of  data  used  in  the  solution.  It  should  be  emphasized  that  this  analysis  is  only 
for  the  single-pass  geometry  with  a  500-naut-mi  near-circuiar  orbit. 

Figure  9  also  shows  the  relative  effect  of  satellite  geometry  upon  navigation  accu¬ 
racy.  Bold  lines  indicate  passes  with  maximum  elevation  angles  of  15°  and  75°  with  12 
to  16  min  of  data  available.  Any  pass  with  maximum  elevation  between  15°  and  75°  and 
7  to  15  min  of  data  has  relative  weight  near  1. 

The  remaining  lines  indicate  only  partial  data  collection,  varying  fror\  minimum 
data  (3  points)  to  half  of  the  pass.  These  relative  weights  may  be  used  to  estimate  a  fix 
accuracy:  (a)  data,  and  (b)  multiplying  by  the  total  system  error,  which  includes  the 
measurement  noise,  satellite  position  uncertainty,  etc. 

The  relative  navigation  error  is  given  as  a  function  of  the  data  span  (or  time)  in 
Fig.  10,  where  it  is  assumed  that  the  maximum  elevation  angle  is  between  15  and  75° 
with  the  data  taken  symmetrically  with  respect  to  the  maximum  elevation  point.  For 
example,  the  figure  shows  a  factor  of  20  in  the  error  distribution  of  a  collection  of  navi¬ 
gations  obtained  using  a  15-min  data  span  versus  a  collection  of  navigations  obtained 
using  a  3-min  data  span. 

It  should  be  noted  that  the  discussion  has  been  limited  to  the  determination  of  three 
navigational  parameters  (latitude,  longitude,  and  a  bias  parameter)  from  a  single-pass 
configuration  which  allows  up  to  16  min  of  data  collection. 


CONCLUSIONS 

It  is  concluded  that,  using  the  techniques  outlined  in  this  report,  (a)  a  navigator  can 
determine  his  position  and  synchronize  his  clock  by  using  the  passive  ranging  technique, 
(b)  the  measured  range  data  may  be  used  to  obtain  a  fix  regardless  of  whether  the  user's 
frequency  is  known  or  unknown,  and  (c)  the  time  difference  may  be  obtained  in  either 
case.  The  doppler  data  can  be  used  to  obtain  a  fix  provided  that  clock  synchronization 
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Fig.  9  -  Single-pass  geometry  associated  with  a 
near-circular  500-naut-mi  orbit.  The  x's  on  each 
trajectory  represent  times  at  which  navigational 
data  might  be  obtained  using  the  satellite.  The 
numbers  are  relative  GDOP  (Geometric  Dilution 
of  Precision)  factors  which  indicate  the  effect  of 
partial  data  collection  during  a  satellite  pass. 


Fig.  10  -  Relative  navigation  error  as  a  function  of 
the  time  during  which  data  is  obtained  (called  the 
data  span).  For  example,  the  relative  error  is  20 
times  as  great  for  a  3-min  data  span  as  compared 
to  a  15- min  data  span. 
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is  obtained  from  an  independent  source,  and  that  clock  synchronization  cannot  be  obtained 
from  integrated  doppler  data.  The  range  and  doppler  data  may  be  used  to  obtain  an  instan 
taneous  fix,  assuming  the  use  of  suitable  techniques.  It  is  further  shown  that  for  the  500- 
naut-mi  single-pass  configuration  the  maximum  elevation  angle,  the  amount  of  data,  and 
the  symmetrical  character  of  the  data  are  important  factors  in  navigation  fix  accuracy. 
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